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We propose a scheme that allows to laser cool trapped 
atoms to the ground state of a one-dimensional confining po- 
tential. The scheme is based on the creation of a dark state 
by designing the laser profile, so that the hottest atoms are 
coherently pumped to another internal level, and then re- 
pumped back. The scheme works beyond the Lamb-Dicke 
limit. We present results of a full quantum treatment for a 
one- dimensional model. 



I. INTRODUCTION 

The ultimate goal of laser cooling of neutral atoms 
stored in a trapping potential is the cooling to the vi- 
brational ground state. This is of interest in the context 
of observing effects related to quantum statistical prop- 
erties of atoms tj and object of present experimental 
investigations J^,3|- The hope is that these efforts might 
ultimately open a second route to Bose-Einstein conden- 
sation of dilute gases 

Laser cooling to the ground state of a trap has been 
achieved experimentally for single trapped ions, us- 
ing sideband cooling techniques || in the Lamb-Dicke 
limit Sideband cooling is based on the selective 

laser excitation of the low-frequency sideband in a har- 
monic trap, leading to optical pumping into the vibra- 
tional ground state. The requirements are that the mo- 
tional sidebands are spectroscopically resolved (which is 
achieved in the strong confinement limit), and that the 
spatial dimension of the ground state oq is much smaller 
the the wavelength of the cooling laser A (in order to 
avoid the heating produced by spontaneous emission). 
This latter condition is represented by the relation 77 -C 1, 
where 77 is the Lamb-Dicke parameter 77 = 27rao/A. In 
one of our recent publications we have discussed 
extensions of sideband cooling beyond the Lamb-Dicke 
regime. Over the last few years other laser cooling tech- 
niques have been developed, which achieve ground state 
cooling in trapping potential and which were originally 
developed for free atoms [pPJioj]. In particular, for free 
atoms they have allowed to achieve temperatures below 
the recoil limit E < Er, corresponding to an atomic 
de Broglie wavelength larger than the wavelength of the 
light A. This has been obtained experimentally, by opti- 
cal pumping into a velocity selective dark state , or by 
cooling with a sequence of shaped Raman pulses where 



the frequency spectrum of the light is tailored so that 
atoms with near zero velocity are no longer excited [ p"2[ . 
These subrecoil cooling techniques as well as sideband 
cooling are all versions of "dark state cooling" . The fun- 
damental idea is the decoupling of a quantum state from 
radiation, and the accumulation of atoms in this state 
by spontaneous emission. The dark state condition can 
be created with various mechanisms; in sideband cooling, 
the trap ground state is dark because it is off-resonant 
from any other state, since the laser is tuned to the red 
motional sideband. The free particle subrecoil cooling 
schemes, instead, use dark states in momentum space. 

In the present paper we will discuss a cooling scheme 
which is based on the creation of a dark state in position 
space with the help of an appropriate spatial profile of 
the cooling laser, so that this state (typically the ground 
state of the trap) is not or only weakly excited by the laser 
p3| . This condition can be realized, for example, using 
the laser in the doughnut mode Jl4| with the axis aligned 
with one of the trap axis. In particular, this cooling 
mechanism works for Lamb-Dicke parameters 77 > 1 , and 
outside of the strong confinement regime (i.e. when the 
sidebands are not spectroscopically resolved). We will 
show that the cooling scheme is quite efficient, allowing 
to cool a significant fraction of atoms to the ground state 
of a trapping potential. The model we will consider is 
one-dimensional, but the scheme is readily extended to 
two and three-dimensional situations. 

The paper is organized as follows. In Sec. IIA, we 
describe qualitatively the cooling scheme. In Sec. IIB 
we develop the Master Equation describing the system 
dynamics, and in IIC we discuss the existence of a sta- 
tionary solution. In Sec. Ill we present the numerical 
results, obtained by taking a laser intensity profile cor- 
responding to a doughnut mode with: (i) atoms trapped 
by an harmonic oscillator, and (ii) atoms trapped by an- 
other confining potential. In this last case, the efficiency 
of the cooling mechanism is improved. 



II. THE MODEL 

A. The laser cooling scheme 

We consider laser cooling in a three-level A-system, 
as illustrated in Fig. 1. We denote by \g) and |e) atomic 
ground (or metastable) states which are connected to the 
excited state |r) by dipole transitions. The atoms are 
confined in a one-dimensional harmonic potential with 
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oscillation frequency v. The quantized trap levels will be 
denoted by |n). 

The laser cooling scheme consists of the repetition of a 
sequence of three pulses. In a first step, atoms in \g) are 
coherently transferred to the state |e) by an off-resonant 
Raman pulse (Fig. la). The purpose of this pulse is to 
excite atoms which are not in the ground state of the 
trapping potential. This is achieved by assuming a spa- 
tial distribution of the laser which has essentially zero 
intensity in a the region where the ground state wave- 
function \g, 0) takes on appreciable values. On the other 
hand, the spatial distribution of atoms in excited trap 
states overlaps with the laser light so that these atoms 
will be transferred to the state |e). The second step con- 
sists of laser cooling of atoms in state |e), thermalizing 
the atomic center of mass distribution to a few recoil 
energies Er (Fig. lb). This is achieved by coupling |e) 
to other internal atomic levels employing one of the stan- 
dard laser cooling schemes (e.g. polarization gradient 
cooling) [^5|. This is used to provide a mechanism for 
confinement of atoms in a finite range of trapping levels,. 
In a third step atoms are optically pumped into \g) using 
a laser tuned on resonance to the transition |e) — > \r) 
followed by spontaneous transition to \g) (Fig lc). We 
assume that the branching ratio between the decay chan- 
nels |r) — > \g) and \r) — ► |e) is large, so that we can safely 
neglect the latter decay channel. Atoms decaying from 
|r) to the ground state \g) will be distributed over a range 
of trap levels of the order of a recoil energy, including the 
ground state \g,0). Repetition of this cooling cycle will 
accumulate atoms in the spatial dark state \g, 0). 

B. Master Equation 

We give here the master equation describing the evo- 
lution of the system during the sequence of pulses. Let 
p be the density matrix describing the system of atoms 
interacting with the laser field. The master equation can 
be written in the form [flG| 

j t p(t) = Lp(t), (1) 

where L is a linear superoperator whose form depends 
on the particular cooling step. The formal solution is 
p(t) = e i '* _t °- ) p(to)- The scheme consists of three steps, 
and the density matrix at the end of a sequence of pulses 
will be 

p(t 3 ) = e^^e^-^e^-^pito), (2) 

where Li is the linear superoperator corresponding to the 
i-th pulse of the sequence, and p(t ) is the density matrix 
at the beginning of the sequence of pulses. We derive now 
the explicit form of the evolution operators. 

First step: we assume that the evolution is coherent. 
We eliminate the excited state |r) in perturbation theory 
and reduce the dynamics to an effective two-level system 



I g), \e). Assuming that the lasers are tuned to exact 
Raman resonance and working in the frame rotating at 
the laser frequency, we write 

= e^ H ^^p(t )e iH ^- t0 \ (3) 

where H is the Hamiltonian describing the evolution dur- 
ing the pulse, 

H = Ho + H h (4) 

with 

H a = fiva)a, (5) 

and 

Hi= MW^ e ik& + &e -ik&y (6) 

Here, a and are creation and annihilation operators 
for the harmonic oscillator, v is the trap frequency, cr — 
|e)(g| and a — \g)(e\ are the dipole raising and lowering 
operators, respectively, fl (x) is the Rabi frequency and 
k the wavevector of the two-photon transition along the 
(cooling) x axis. Denoting by At = t\ — to the time 
duration of this first pulse, we will assume that is At <C 1, 
so that we can safely neglect Hq in (|J) 

p{t x ) « c-^ At p(t a )e^ At . (7) 

Second step: Laser cooling of atoms in state |e) via 
an auxiliar level [see Fig. lb] destroys the coherences be- 
tween |e) and \g), while atoms left in \g) evolve according 
to (||) . We do not specify the cooling mechanism, assum- 
ing only that at the end of the pulse the atoms in |e) are 
thermalized, and described by the density matrix 

\e)(e\®p* = \e){e\® ^f^^ , (8) 

where N is related to the average number of vibrational 
quanta after this laser cooling process. Using the prop- 
erties of (7^ , a in (^), we have 

p(t 2 ) = e-^° T - \g) ( g \ cos[n(x)At/2]p(t ) (9) 
x cos[n(f)At/2] \g)(g\e m ° T ^ +C|e)p th (e|, 

where T scp = t 2 — t\ is the duration of this step, and £ is 
the probability of occupation of |e), i.e. 

C = trace{|e)(e|p(i 1 )} (10) 
= 1 - trace{|g) (g\ cos [Q.(x)At/2] p(t ) cos [fi(x)At/2]}. 

Third step: the thermalized atoms in |e) are optically 
pumped into \g). The linear superoperator describing 
the master equation is 

p 

L3P = -i[H , p] - — (pa '<7 + a <ap) (11) 
+ T f duN{u)ae lk ^ 1+u)i ! pth e -i*,(l+u)«-t . 
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with r the optical pumping rate, N(u) the angular dis- 
tribution of the emitted photons, and k p wavevector 
of the pumping laser that propagates along x, k p = 
2ir/\ = r) /ao- Assuming u ^> £3 — t%, we can ne- 
glect the free evolution during this pulse. Furthermore, 
for £3 — ti ^> r _1 all the atoms in |e) are pumped in \g) 
during the pulse. Taking that into account, we find: 



p(t 3 ) = e- liloi ^\g){g\cos[Q(x)At/2}p(t ) 



(12) 



x cos[n(x)At/2]\g)(g\e iHoT ^ 
+ Cj duN{u)e lk ^ 1+u ^\g)p th (g\e~ lk ^ 1+u ^. 

We can now project onto \g), and define p g = (g\p\g). De- 
noting by F = j\ duN(u)e lk " il+u ^p th e- lk '' (1+u ^ the 
matrix describing the feeding contribution, we finally 
have 



p g (t 3 ) = e' lHoT ^ cos[n(x)At/2]p g (t ) 



(13) 



x cos[fl(x)At/2]e 



iH T 3af 



(F 



This gives a mapping between the density operator at 
the beginning and at the end of a given sequence of 3 
pulses. From the above expressions, we can deduce that 
for large traps (i.e. for large values of the Lamb-Dike 
parameter) the cooling will slow down and become less 
efficient. The reason for that is twofold: on one hand, 
the kick provided by the optical pumping process will 
distribute the atoms among a wider range of states. On 
the other hand, since we have assumed that the thermal 
distribution p th has a mean energy of the order of the 
recoil, the mean occupation number will be larger ((n) sa 
E R /hv = V 2 ). 



C. Stationary state 

Ideally, we would like to use the cooling mechanism to 
cool all the atoms to the ground state (or any other pure 
state). Therefore, we have to analyze the shape of il(x) 
which gives rise to the optimal cooling. In this subsection 
we analyze qualitatively the conditions for the spatial 
profile of the laser £l(x) to provide an efficient cooling. 

We then look for stationary solutions to the map- 
ping (13) which corresponds to pure states. A neces- 
sary condition for that is C = 0. Then, we substitute 
C = in ([[3]) and impose the stationary regime, i.e. 
Pg{h) = Pg{to) = Pg, with p g = |/)(/|. We obtain that 
these conditions are satisfied if the following equation al- 
lows for a solution: 



x c 



iH T s < 



\f)e i< f> = cos[n(x)At/2](x\f), 



(14) 



with <f> arbitrary phase. Let us denote by R the re- 
gion where (x\f) is different from zero. It can be easily 
shown that (Q) implies that in the region R the rela- 
tion I cos[f2(x) At/2] I = 1 must be fulfilled. Physically, 



this means that if we want to have a dark state |/) its 
wavefunction must be completely distributed in the spa- 
tial region R where there is no laser excitation (that is, 
where the atoms perform complete Rabi oscillations) . On 
the other hand, we want to have only one dark state, and 
therefore other states must have a finite excitation prob- 
ability. This means that they must be nonzero outside 
R. Therefore, for achieving cooling to the ground state, 
one should design the intensity profile of the laser in such 
a way that only in the region where the ground state is 
localized, the laser accomplishes no excitation. This con- 
dition can be only approximately fulfilled in a harmonic 
trap, where the ratio among the spatial dimensions of 
the ground state and the first excited state is 1/V2- In 
this case, we expect to not be able to cool all the atoms 
into the ground state. Cooling into the ground state will 
be possible in the case of atoms trapped by a confining 
potential, which better localizes the ground state respect 
to the first excited state. 

The qualitative argument given above does not de- 
pend on the separation time T sep in the second pulse. 
In reality, the selection of this time can play an im- 
portant role, especially in the case of a harmonic po- 
tential. To see that, let us consider the particular case 
vT scp = 21tt, with I integer. Then, according to ( |l3| ) the 
positions x in the region R in which | cos[f2(a;) At/2] | 7^ 1 
will be emptied, whereas the atoms will be accumu- 
lated in the regions where | cos[f2(ir)Ai/2]| = 1. In fact, 
from (14) we see that the state |/) is a stable solu- 
tion if the conditions cos(Vt(x)At/2)\f) = exp(i0i)|/) 
and if exp(— iHoT scp )\f) = exp(i02)|/) are fulfilled, with 
01+02 = 4>- For vT scp = 21-k this means that each 
wavepacket completely distributed inside R at t is so- 
lution (i.e. a dark state), since in T sep it undergoes a 
full oscillation and therefore recovers the original form. 
In that case, laser cooling will not be possible since there 
will exist many wavepackets that will remain dark. In the 
case vT sep = (21 + 1)tt, one can consider two subspaces 
St ={\0), |2), |4), ...} and S 2 ={\1), |3>, |5), ...}, in which 
the evolution during the second pulse is periodic. That 
is, if we can form a wavepacket as a linear combination of 
the states of Si or of S2 which is spatially distributed in 
R at to, this wavepacket will be a dark state. Hence, laser 
cooling will not be possible either. The argument can be 
repeated now for times fulfilling vT scp — mr/m, and sub- 
spaces can be constructed so that wavepackets belonging 
to them will be dark states, and therefore laser cooling 
will not be possible. We can overcome this problem very 
easily by choosing random separation times T sep so that 
the only wavepacket that is a dark state is precisely the 
ground state. 



III. RESULTS OF CALCULATIONS 

In the following, we show results obtained assuming 
a doughnut mode as spatial distribution of the Raman 
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pulse. The reason for this choice is that it provides a flat 
profile in a finite region which, as we have shown above, 
is required to obtain a suitable cooling. In particular, 
we discuss the parameter regime of the laser for which 
the cooling efficiency is optimal. Finally, we explore the 
role of T scp in the cooling dynamics. We analyze first the 
case of a harmonic trapping potential and then another 
potential which is best suited for laser cooling using the 
present scheme. 

We take as spatial distribution for the Raman pulse 
laser the following class of functions: 



/(a 



2d 



exp 



x/oq 



(15) 



Here a determines the spatial width of the pulse in units 
of the harmonic oscillator ground state width ao, and the 
exponent n is an integer number, which denotes the order 
of the doughnut mode and thus determines the flatness 
around x = 0. In the following, we explore the depen- 
dence of the cooling efficiency on the parameters n and 
a, for atoms whose center of mass is trapped by a har- 
monic potential with frequency v. We take the Rabi fre- 
quency Q(x) oc /(£)■ Furthermore, the projection along 
the cooling axis (x) of the two-photon wavevector appear- 
ing in (||) is k = 2n/X, where A is the wavelength of the 
transition \r) — > |e), since the doughnut mode photons 
propagate in a direction orthogonal to the x axis. 

In Fig. 2a and Fig. 2b we show the population P° = 
{5, 0|yO g |<7, 0) after 1500 (curve with o) and 2500 (curve 
with x ) sequences of pulses as a function of the width a 
and of the exponent 2n, respectively. We have obtained 
these results using Eq. (13), with a Lamb-Dicke parame- 
ter rj = 2TTagJX = 5, and with a thermalized distribution 
p th in Eq. (g) with N = rj 2 = 25. From Fig. 2a the exis- 
tence of an optimal width a for a given exponent (in the 
present case: 2n = 4) is evident. The optimal a depends 
on n, since it must fulfill the requirements of Eq. (|14|). In 
Fig. 2b we plot P® as a function of 2n, where the values 
of the width have been optimized. We see that for n > 4, 
around 80% of the atoms are cooled to the ground state. 
Figure 3a shows the population P™ = (n\p g \n) as a func- 
tion of the vibrational number n after 2500 sequences of 
pulses for n = 4 and a — 4, and Fig. 3b shows P® as 
a function of the number of pulses. In the inset of this 
figure the final spatial distribution is shown (solid line), 
while the spatial distribution cos[fi(a;)Ai/2] is plotted as 
a dashed line for comparison. In all the cases, the final 
density matrix is not a pure state but a mixture, diagonal 
in the number states basis. According to the discussion in 
the previous Section, we have chosen random separation 
times. We have also checked that by fixing the separation 
times T sep , the cooling efficiency decreases significatively. 

The efficiency of the cooling scheme in populating the 
ground state depends on the difference between the reex- 
citation rate of the vibrational ground state and the ones 
of the first few excited states. To the extent that one 
can design a confining potential so that the ground state 



is much more localized in space than the other eigen- 
states, the cooling efficiency will improve correspond- 
ingly. A confining potential satisfying the above require- 
ment could be experimentally realized, for example, by 
adding to an harmonic trap an additional optical poten- 
tial. In the following we will assume a potential 



V{x) = -2mu 2 



1 + gx 2 



(16) 



where e and g characterize the shape of the potential 
near the center. An example of how the potential ( fl6| ) 
looks like is shown in Fig. 4a. Figure 4b shows the spa- 
tial widths of the first twelve eigenstates for g = 2000/oq 
and e — 19 x 10 5 (curve with o), compared with the 
corresponding widths of the harmonic oscillator eigen- 
states with the same v (curve with x ) . For this choice of 
the parameters e and g, the eigenvectors of the potential 
tend to the number states of the harmonic oscillator for 
n > 1. In the present example, the ground state has spa- 
tial width ao/5, and it is thus well localized respect to the 
other eigenstates. The laser cooling scheme works very 
efficiently in populating the ground state, as it may be 
seen from Fig 5: here the population of the ground state 
in the trapping potential V(x) is plotted as a function 
of the number of pulses (solid line) for a doughnut mode 
with a = 8 and 2n = 2 [c.f. Eq. @]. The Lamb-Dicke 
parameter for the eigenstates n > 1 and for the internal 
excited state |e), is r] — 5, while the ground state has 
spatial width corresponding to an effective Lamb-Dicke 
parameter rj = 77/5 = 1. We can see that during the 
cooling almost all atoms are accumulated in the ground 
state of V{x). The dashed line in Fig. 5 shows for refer- 
ence the evolution of the population of the ground state 
of a pure harmonic oscillator potential with r\ = 1. In 
this second case, we have used a doughnut mode with 
2n = 2 and a — 10. Note that the sample with 77 = 1 
is colder at the beginning, since its initial distribution is 
thermalized with (n) « 4. 

The result obtained here is rather general and can be 
summarized as follows: through the creation of a more lo- 
calized ground state, obtained by a suitable perturbation 
of an harmonic potential of a given 77, the present cool- 
ing mechanism allows to cool all the atoms to the ground 
state, even for the Raman pulse intensity distribution in 
a low order doughnut mode. 



IV. CONCLUSIONS 

In this paper we have studied laser cooling of trapped 
atoms to a single quantum state, based on creating a 
dark state in position space by reexciting the atoms with 
cooling laser which has no spatial overlap with the trap 
ground state. The scheme can be applied in the limit 
of weak confinement (r > v) and outside the Lamb- 
Dicke regime. Assuming a doughnut mode as a model 
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of the spatial profile of the laser, our numerical calcula- 
tions have shown that more than 80% of the atoms can 
be cooled to the ground state of a harmonic oscillator po- 
tential. Furthermore, this number can be increased sig- 
nificantly for trapping potentials where the ground state 
is much better localized than the first few excited vibra- 
tional states. The randomization of the separation time 
in the harmonic oscillator case is fundamental in order 
to achieve laser cooling to the ground state. 

The model presented here is one dimensional, but it 
can be easily extended to two dimensions. In that case, 
one could use the plane orthogonal to the doughnut mode 
axis obtaining similar results as in the one dimensional 
case. The reason for that is that the heating processes 
(due to photon scattering during repumping) that we 
have considered in the one dimensional case will be simi- 
lar to the two dimensional one. Experimentally, the one- 
dimensional and two-dimensional schemes correspond to 
the employement of a single doughnut mode, whose axis 
coincides with the one and two-dimensional trap axes 
respectively, and which dark region is shaped into the 
size of the one-dimensional and two-dimensional trap 
ground state respectively. The extension to three dimen- 
sions could be achieved crossing two doughnut modes, 
which propagate in directions that are orthogonal one to 
the other, and whose axes cross at the center of a three 
dimensional trap. In this way a three dimensional dark 
region would be created, that could be shaped into the 
size of the trap ground state. The scheme may be useful 
for achieving Bose-Einstein condensation with laser cool- 
ing, and more in general in preparing with laser cooling 
a macroscopic occupation of a state of motion. 
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FIG. 1. Internal configuration of the atom, and scheme 
of the sequence of pulses. Internal configuration: \g), \e) sta- 
ble or metastable states, \r) excited state, with \g) — ► \r), 
\e) — > \r) dipole transitions. Scheme of the sequence: (a) co- 
herent Raman pulse, which corresponds to a doughnut mode 
on \g) — > jr); (b) cooling to the recoil on the level |e), obtained 
coupling \e) to a fourth internal level, (c) optical pumping 
into \g), obtained tuning a laser on resonance on |e) — ► \r). 
We assume that the branching ratio between r) — > \g) and 
\r) — ► |e) is very large, so that the decay along \r) — > \e) can 
be safely neglected. 

FIG. 2. Plots of the efficiency after 1500 (curve with 
o) and after 2500 (curve with x) sequences of pulses, in 
function of (a) the width a, for 2n = 4 and (b) the ex- 
ponents 2n — 2, 4, 6, 8, with optimized widths, respectively 
a = 30,4.2,3,2.2. In all cases max(fi(a;) At/2) = 0.6tt. uT scp 
is random, taken from a flat distribution that varies in the 
interval [0.1, 1.1]. 

FIG. 3. (a) Population in function of the vibrational num- 
ber state after 2500 sequences of pulses, for a harmonic oscilla- 
tor with ij = 5 and Raman laser corresponding to a doughnut 
mode with 2n = 4, a = 4 and max(f2(x)At/2) = 0.67T. In- 
lay: spatial distribution of the sample after 2500 pulses (solid 
line) and cos(S7(x)At/2) (dashed line) plotted for compari- 
son, (b): Ground state of the trap population in function 
of the number of pulses, under the same condition as above. 
The initial distribution and the excited state distribution af- 
ter the recoil cooling are described by p th with N — 25. vT scp 
is random, taken from a flat distribution that varies in the 
interval [0.1, 1.1]. 



FIG. 4. (a): plot of V(x) for arbitrary values of e and 
g. (b): (Ace 2 ) 1 '' 2 in function of the first 12 eigenstates of 
V(x) (curve with o) and of the first 12 ones of the Harmonic 
oscillator potential ( curve with x), for e = 19 x 10 5 and 
g = 2000/a 2 ,. 

FIG. 5. Ground state occupation probability versus the 
pulse number. Solid line: ground state of atoms trapped by 
V(x), with rj = 5 and rj' — 1, p th with N = 25 and Raman 
laser corresponding to a doughnut mode with 2n — 2 and 
a = 8. Dashed line: ground state of an harmonic oscillator 
with ij = 1, p th with N = 4 and Raman laser corresponding 
to a doughnut mode with 2n = 2 and a = 10. In both cases 
max(f2(a;)At/2) = 0.67T. uT scp is random, taken from a flat 
distribution that varies in the interval [0.1, 1.1]. 
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